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Abstract 

The aim of this paper is to introduce an insurance model allowing reinsurance and 
dividend payment. Our model deals with several homogeneous contracts and takes 
into account the legislation regarding the provisions to be justified by the insurance 
companies. This translates into some restriction on the (maximal) number of con- 
tracts the company is allowed to cover. We deal with a controlled jump process in 
which one has free choice of retention level and dividend amount. The value func- 
tion is given as the maximized expected discounted dividends. We prove that this 
value function is a viscosity solution of some first-order Hamilton-Jacobi-Bellman 
variational inequality. Moreover, a uniqueness result is provided. 
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1 Preliminaries 



A common problem of the insurance companies is to find a strategy allowing 
to satisfy the claims appearing either from the insured parties as consequence 
to specified peril or from the shareholders in terms of dividends. To reduce 
their risks and protect themselves from very large losses, the companies usually 
choose to pay some of the premiums to a third party. This process is called 
reinsurance, and it commits the third party (the reinsurance company) to 
cover a certain part of the claims. It is obvious that the insurance company 
controls the contracts to be reinsured as well as the dividends to be paid to 
the shareholders. These elements justify the framework of stochastic control. 
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This paper considers a utility function given as the maximized expected dis- 
counted dividends. In the literature, this approach has been first used by 
Jeanblanc, and Shiryaev (1995). In their model, the capital of an insurance 
company is described with the help of a standard Brownian motion and the 
dividend payment strategy is understood as control process. More precisely, 
they deal with the following model 



where \i and a are arbitrary constants, W is a 1-dimensional standard Brow- 
nian motion and Z is an adapted, non decreasing, right-continuous process 
which represents the dividend payment strategy. 

In Asmussen et al. (2000), a model concerning excess-of-loss reinsurance and 
dividend payment has been studied. They use diffusion and proportional rein- 
surance for their model. More exactly, they take as model for the capital of 
the insurance company the process given by the following equation 



where < a t < 1 stands for the retention level. In the case where the rate of 
dividend pay-out is unrestricted, they characterize the value function as the 
(classical) solution of some associated Hamilton- Jacobi-Bellman equation. 

The same problem is studied by Mnif, Sulem (2005), but the claims are rep- 
resented by a compound Poisson process. In their collective risk model, a 
retention level is an adapted process at which specifies that, for a claim y, the 
direct insurer covers yAat, while the reinsurance company covers the remain- 
ing (y — a t ) + . They consider a single insurance contract and the reserve of 
the insurance company satisfies 



where \i is the random measure associated to the compound Poisson process. 
In the above equation, p(a t ) is the actual premium of the insurance company 
given the retention level a. The process L describes the pay-out of dividends for 
shareholders and it is an adapted, cadlag process such that L t — L t - < X t - for 
all t > 0. The value function is defined as the maximized expected discounted 
dividends until the ruin time r, 



here r is some positive discount factor. The authors proved that, under the as- 
sumption that the value function satisfies the dynamic programming principle, 
V is a viscosity solution of the associated Hamilton- Jacobi-Bellman variational 
inequality. 



dX t = iidt + adW t — dZ t 



dX t = a t (jidt + crdW t ) - dZ t 



dX t = p(a t )dt — (y A a t ) fi(dtdy) — dL, 



B 



V{x) = sup E 

(u,L) 
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In the present paper we consider the problem of optimal reinsurance and 
dividend pay-out with several insurance contracts. We will prove that in the 
framework of the collective risk model, even if the invested initial capital 
is arbitrarily small, one can expect a gain which exceeds an a priori fixed 
positive constant. Indeed, this comes from the fact that, independently of its 
initial capital, the model allows the insurance company to sell one contract. 
However, as it is precised in section 2, in the case of insurance companies, the 
codes of law impose that, at any time, these companies should be able to justify 
enough resources to cover the obligations contracted towards their clients. This 
condition imposes an upper limit for the number of contracts the company 
can have. In the work we present here, several contracts are considered. We 
obtain a stochastic differential equation with respect to a random measure 
and introduce the utility for the shareholders as in Mnif, Sulem (2005) to be 
the maximized discounted flow of dividends. We prove that the value function 
is regular enough (enjoys the Lipschitz property) and satisfies the associated 
Hamilton- Jacobi-Bellman Variational Inequality in the viscosity sense. We also 
provide an uniqueness result for the viscosity solution in the class of continuous 
functions of at most linear growth. We emphasize that the limitation of the 
number of contracts which comes from the codes of insurance, allows us to get 
the Lipschitz property of the value function V. This property insures that an 
initial capital close to will induce a zero-expected gain (unlike the collective 
risk model). Moreover, in this case, the dynamic programming principle follows 
in a standard way, while it was only assumed by the authors of [9]. 

The paper is organized as follows. In the first section we present a simple 
example showing the limits of the collective risk model. The second section 
is concerned with the insurance problem with several contracts. We introduce 
the model, the basic assumptions and prove some elementary properties of 
the value function V. In the third section, we show that the value function 
is a viscosity solution of the associated Hamilton- Jacobi-Bellman variational 
inequality. The fourth section provides a comparison result which allows to 
obtain the uniqueness of the viscosity solution for the given variational in- 
equality. A numerical example is given in the last section. 



2 The limits of the collective risk model. A counter example 

We consider the following special case of the collective risk model introduced 
by Mnif, Sulem (2005). We assume that the claims are generated by a Poisson 
process N with intensity 1 on a complete probability space (f2,jF, P). We 
denote by (Ft) t>Q the filtration generated by the random measure associated 
to N, completed by the family of P-null sets. Given an T t — adapted process 
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a t G [0, 1] (retention level), the premium rate is 

p(a t ) — ki — k 2 + (1 + fo) a t , for all t > 0, 

where < k\ < k 2 are proportional factors. Moreover, if L denotes the F t - 
adapted process of cumulative dividends, then the reserve of the insurance 
company satisfies the equation 



The process L should be right-continuous, non-decreasing and such that L - — 
and L t - U_ < X^ L for all t > 0. We introduce the first jump time for 
the Poisson process iV 



Obviously, t± is of exponential law with intensity 1, and, in particular, 




7i = inf {t > : N t = 1} . 



P(ri>l) 



-i 



If we consider the strategy (a, L) given by 



a 



= 1, 



k L t (w) = J{ ri >i}(^)/{t>i}(t), 



then (a, L) is admissible and the ruin time 



> 1 on {n > 1} . 



Indeed, 




and on {j\ > 1} we have that 



x + {\ + kx)t, 



for all t < 1. 



It follows that 




for all a; > 0. Obviously 



v (o+) > e -( r+1 ) > 0. 
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Therefore, investing an arbitrarily small capital in the insurance company, we 
expect to gain more than e~( r+1 '. This contradicts theorem 3.3 in [9]. This 
problem is due mainly to the fact that, independent of the initial capital, the 
insurance company is allowed to hold one contract. 

However, the insurance law requires that, at any moment, the companies 
should be able to cover any liabilities that have been incurred on insur- 
ance contracts as far as can be reasonably foreseen. Experience of similar 
claim development trends is of particular relevance. Usually, the solvency 
margin is computed with respect to both the premium rates and the aver- 
age claim. According to the current Solvency I prudence regime, "the life 
insurance capital requirements are arrived at by multiplying a factor of 4% 
to the mathematical reserves of participating business (for unit-linked busi- 
ness the factor is reduced to 1%) plus a factor of 0.3% to the sum- at-risk" 
(CEA and Mercer Oliver Wyman, Solvency Assessment Models Compared, 
http: / / www.cea.assur.org/ cea/ download/ publ/ article221.pdf [ ) . 

The suitable formulae should take into account the specificities of life, non-life 
and reinsurance business. Various methods are, therefore, available. To give 
an example, according to the French legislation (Code des Assurances, R334- 
13) for the life insurance, the solvency margin (to be replaced by the Solvency 
Capital Requirement for Solvency II) should be superior to the result obtained 
by multiplying 0,3% of the capital under risk with the ratio between the capital 
under risk after reinsurance and the capital under risk before reinsurance 
computed for the previous exercise. The latter ratio cannot be inferior to 50%. 
To keep it simple, at time t the result obtained by multiplying a constant Co 
(depending on previous experience and the type of insurance business) by 
the average claim per contract and by the number of contracts n t should not 
exceed the fortune of the insurance company: 

(o x fi ( x average claim < fortune at time t. (1) 

Corroborating these elements, it appears obvious that the simple collective risk 
model should be improved to a model involving several contracts. We empha- 
size the fact that only quantitative requirements are taken into consideration 
(therefore, the model covers only part of Solvency II Pillar 1 requirements). 



3 The insurance problem with several contracts 

We introduce a complete probability space (ft, J 7 , P). In order to model the 
claims, as for Mnif, Sulem (2005), we use a compound Poisson process given 
by a random measure fi(dtdy) on IR + x B, with B C K + \ {0} . Moreover, 
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we assume that the compensator of p takes the form dt7r(dy) and that the 
measure 7r is finite 7r(dy) = (3G(dy) for some probability measure G(dy) on B 
and some positive constant j3. 

Throughout the section, we let Y denote a generic random variable distributed 
according to G(dy). 

We consider the natural filtration (J 7 t)t>o generated by the random measure 
/i. We call retention level any (jF 4 )-adapted process (u t ) t > which specifies 
that, given a claim y at time t > 0, the direct insurer covers y A u t while the 
reinsurance company covers the excess of loss (y — u t ) + . 

Since we are going to consider several insurance contracts, we introduce a 
function / depending both on the number of insurance contracts and on the 
risk taken by the company to model the claims / : R + x R — ► M. + . If the 
company chooses some retention level u t , then the actual premium rate per 
contract is given as in Asmussen et al. (2000), or, again, in Mnif, Sulem (2005) 

p(ut) = (1 + h)(3u - (1 + k 2 )(3E [/(l, (Y - u t ) + )} for all t > 0, (2) 
where fej are real constants satisfying < k\ < k 2 and 

v=l f(l,y)G(dy) = E[f(l,Y)\. (3) 

JB 

The first term in is the premium received from the client, while the second 
term is the quantity paid to the reinsurer. 

Given the initial fortune x > and the retention level u, if L stands for the 
{Ft)— adapted process representing the cumulative dividends paid up to the 
time t, n t denotes the number of contracts of the insurance company at time 
t, and Xt ,u ' L the fortune of the company , then we have 

Xt' u,L = x+[ n s p(u s )ds- [ [ f(n s ,y Au s )/J l (dsdy) - [ dL s . (4) 
Jo Jo Jb Jo 

If we denote by a the quantity 

1 



then, from ([T]) we get that the maximum number of insurance contracts is 
H 



n max _ a Xt ,u ' L . We have the following equation 



X^ L = x + a ( X^ L p(u s )ds- I I f(aX^ L ,yAu s )fi(dsdy)- f dL s , 
Jo Jo Jb Jo 

(5) 
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and introduce the cost functional 



/ 

Jo 



•T 



J(x, u,L) — E 



e~ rs dL 



s 



(6) 



where r is some discount factor and r is the ruin time 



T = 



inf {t > : X?' U ' L < 0}. 



Our value function V will be defined as the maximum over some family of 
admissible couples (u, L) of the cost functional J. 

In practice, whenever the solvency condition is not satisfied, one of the follow- 
ing two events may occur. In the first capital infusion from the share- 
holders intervenes. In the second one, an external referee solves the problem: 
either by transferring some of the contracts to other insurance companies, or 
by dissolving the contracts in final phase. The Solvency II framework states 
that as soon as the Solvency Capital Requirement (SCR) is not satisfied, su- 
pervisory action will be triggered. However, if the Minimum Capital Require- 
ment (MCR) is not satisfied, the control authority can invoke severe measures 
(including closure of the company). From the mathematical point of view, we 
do not allow capital infusions, these being obtained by taking a larger initial 
reserve. On the contrary, the latter events may appear and they allow the 
variation of the number of contracts. 

Let us now return to the function / modelling the claims. It is natural to 
suppose that the claims increase with the number of contracts and are null if 
the company has no contract. Moreover, the claims should increase with the 
risks covered and should be if dealing with no risk. If the number of contracts 
is positive and the risk covered by these contracts is not null, then the claims 
are expected to be strictly positive. An utility function is usually supposed to 
be concave. If we are given a concave function v such that v(0) = 0, then 



for any A < 1. Since any nonlinearity in (jSJ) may only come from /, in order 
to obtain the previous property for our utility function V, one should assume 
that / is convex in the first variable. These assumptions give 



v(Xx) > Xv(x), 



Assumption 1 (Al) Suppose that the function f : R + x R 



R + satisfies: 



f{-,y) is convex, non decreasing and f(0,y) = for all y e M + ; 



f(x, •) is increasing and f(x, 0) = 0; 



f{x,y) > if x > and y > 0; 



- f is uniformly continuous on R + x R; 
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- f(x,y) is Lipschitz in x, uniformly in y e M+. 

One expects to cover expenditures through the premium received 

p(u t )> PE[f(hYAu t )\. 

Recall that p(0)-/3E[f(l, 0)] < and that Zinv-oo (p(«) - PE[f(l, Y Au)])> 
(recall the definitions (j2J) and (j3J) of p and z/, respectively) and we obtain 
the existence of some u > such that 

P («)>^[/(i,yA«)], (7) 

for all u > u. Thus, we are going to consider only the retention levels u t 
satisfying 

u t > u- (8) 

One should impose that the dividends paid at some time t do not exceed the 
reserve at the same time. Therefore, we call admissible strategy the couple of 
(J-'t) — adapted processes (u,L) such that u satisfies (JHJ) and L is cadlag, non 
decreasing, L _ = and L t — L t ~ < X^' L for almost every (t, u). We should 
first prove the existence of such admissible strategies. 

Remark 2 If I is an (Ft) — adapted process which is cadlag, non decreasing, 
Zo- = 0, then, for any initial condition x > 0, and any (J-'t) — adapted processes 
u which satisfies there exists a unique Tt~ adapted right- continuous process 
Xt ,a ' 1 with left-hand limits which satisfies the equation 

X?' u > l = x + a I X x s ^ l p(u s )ds- f t+ f f{aX^' l ,yAu s )fi(dsdy)- f dl s (9) 
Jo Jo Jb Jo 

(see also Ikeda, Watanabe (1989) IV, Theorem 9.1). We define the ruin time 
t = inf |t > : Xt' u ' 1 < oj. Obviously, on {t < r} we have Al t = l t — l t _ < 
X^' 1 . Let us define the process 

L t = l t l {t<r} + (Al t AX?l u > 1 ) l {t=T} . 

We get an (J-'t) — adapted process which is cadlag, non decreasing, and Lo- = 0. 
Let X X,U,L denote the solution of |P)) with L instead of I. We notice that (u, L) 
is admissible in the sense that L t — L t - < for almost every (t,u). 

For all initial reserve x > 0, we denote by A(x) the set of admissible strategies 
described above. The value function is defined by 

V(x) = sup J(x,u,L). 

(u,L)GA(x) 

Proposition 3 (Comparison for solutions of ^)) Given two (J-'t) — adapted 
processes u and I such that u satisfies (EJj and I is cadlag, non decreasing, and 
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/o- = 0, and two initial states < x < x', the solutions of (TJ|) X x ' u ' 1 and X x ' u ' 1 
starting from x (respectively x') and associated with the pair (u, I) satisfy 

j£X,u,l -^x ,u,l^ a ^ j- p_ a s 



PROOF. Let us consider the sequence of functions (fi n G C 1 (M) such that 
<f) n {x") = for all x" < 0, < <j/ n (x") < 1, for all x" G R, and <j> n {x") T (x") + 
as n — > oo. A simple application of Ito's formula yields 



(-y-X,U,l v x',U,l\ j | j 

A t — A t J — ii + 1 2 , 



(10) 



where 



a = j\ P {u s ) (x* s ^ 1 - x:'^ 1 ) <p' n (x:^ 1 - x?^) ds, 
h = 

J j 4> n (x x ^ 1 - XsJ u ' 1 - f (aX x ^'\ y A w s ) + / (aX x 2 u '\ «/Au s )) (i(dsdy) 
- J o J b 4>n (x x '^ 1 - X x '_l u ' l ^j jj(dsdy). 



It is obvious that 



h < c J* (x x s ^ l -x : 



x',u,l\ + 



ds, 



where C is a constant independent of x and x'. Since a can be chosen arbitrarily 
small (for that, it is enough to recall a = and then choose an arbitrarily 
small monetary unit such that the quantity v becomes large), we may assume 
that aK < 1 (here K denotes the Lipschitz constant for /). Then the function 
x i — > x — f(ax, y) is increasing for all y G K + . Therefore, we get 

h < o. 

Combining the two estimates for Ii and we have 



E 



< C E 



ds. 



We allow n —>■ 00 to obtain 
E 



f x x, u ,i _ x x',u,iy < C f E f X x,«,i _ X x',u,l 

J «^ L 

Finally, Gronwall's inequality yields 



ds. 



E 



^j£X,u,l X X ' ,U '^ 



0. 



The proof of our Proposition is complete. 
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If the initial fortune is fixed, then the company has to make a choice over 
some family of admissible strategies. One may naturally wonder whether the 
same strategies are valid when dealing with a greater initial reserve or not. 
The answer is affirmative as proven by the following Proposition. 

Proposition 4 // < x < x' are two initial capitals and if (u, L) is an 

admissible strategy for x , then (u,L) is also admissible for x' . 



PROOF. Indeed, if Xf' u ' L (respectively X? ' U ' L ) denote the solutions of 
starting from x (respectively x') associated with the control pair (u,L), then 
the comparison result yields 

x x,u,L < X f,u,L^ dup _ ae Qn ^ ^ x Q 

Now, since L is admissible for x, we have 

L t - L t _ < X?l u ' L < Xf> L , dtdP - a.e. 



and L is again admissible for x' . Moreover, if r denotes the ruin time for X t 
and t' denotes the ruin time for X? ,U,L , then, obviously 

t < t', P — a.s. 



x,u,L 



As one expects, using the previous results, we find that the utility function of 
the insurance company increases with the initial reserve. Since our strategy 
involves a dynamic programming approach, we would like to have finite value 
function. We suppose that the following assumption holds true 

Assumption 5 (A 2) The discount factor r in (0|) satisfies 

2(l + fci)/3 

r > 

Co 

Given an economic framework in which the discount factor r is fixed, the above 
assumption says that the time between two claims is great enough to justify 
the demand for small solvency translated in the small constant fo- 
under this Assumption, we provide an upper bound estimate as well as Lips- 
chitz regularity of the value function. 

Proposition 6 The value function V is non decreasing, enjoys the Lipschitz 
property and satisfies 

V(x)<Kx, (11) 
for some large enough positive constant K. 
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PROOF. The first assertion is straightforward from the previous Proposi- 
tion. In order to establish the upper bound (fTT|) . we notice that 

x ,,u,L <x (1 + *!)/? f'^u,^ 

(o io s 
for all t > 0. Gronwall's inequality yields 

(l+fcl)£ 

X^' L < xe^^\ (12) 

We write Ito's formula for e - rt x x,u ' L and use (fT2|) together with (A2) to 
obtain 

J(x, w, L) < Cx. 

Here C is a constant which may change from line to line. Let us fix x, x' > 0. 
Suppose that (u, L) G A(x + x') and notice that, in this case, (u, -^jLj G 
Indeed, 

x x +x >,u,L = ( x + x ^ + a r* X X s +X '> U > L V (u s )ds 

Jo 

- / f(aX x + x ' u,L ,y Au s )fi(dsdy) - dL S} 
Jo Jb Jo 

and, by multiplying the latter equality by ^tp, we get 

' x x+x '' u ' L ■>■ I f a,,{u..)-!— f .X ds 
Jo x + x' 



x + x' 



- I I J ( aX s- X ' ,U,L , y A u s) ^{dsdy) 
Jo Jb x + x' K ' 



* 1 x 



d L 



o \x + x' 



On the other hand, 



x,u, x , L f l , . x.u, x , L 

X t x+x =x+ ap{u s )X s x+x ds 
Jo 

-jf o J B f(aX.- ,yAu s jfi(dsdy) 



t ' x 



d[ L. 

o \x + x' 

Now, let the functions <f> n G C 1 (R) be such that 4> n (x") = for all x" < 0, and 
< <f)' n {x") < 1, for all x" G R, and <f) n {x") | {x") + as n -> oo. We make the 
following notation 

— ,u,L y 

x + x' 

We apply Ito's formula to have 

Y t -X*' U '^ L ) =h + I 2 , (13) 
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where 



h = J ap(u s ) [Y s -X. ^ J <j>' n [Y s - X s ) ds, 

h = 

J b K (Y s _ - XZ^ L - ^/ yAu s )+f (aX X s L U ^\ t/A« s )) ^dsdy) 

- ^ j B (f) n [Y s _ - X s _ =^ j //(dsdy). 

It is obvious that 

h<C [Y s -Xs OT ds, 



jo \ / 
where C is a constant independent of x and x' , and we use the convexity of / 
in the first variable and /(0, ■) = 0, together with the monotonicity of (f) n to 
get (as in the proof of the comparison result), 

h < 0. 

Thus we obtain, as in the comparison result, 

■'■ x x+x>,u,L < X -^7W L dup _ a ^ on [Q) oo) x n 



X + x' 



Obviously, (u, -^jLj is an admissible strategy for the initial reserve x. If r 
is the ruin time for the strategy (u, L) for the initial reserve x + x', then the 
above inequality states that the ruin time for the strategy (u, -^jL^j when 
the initial reserve is x is greater than or equal to r. Therefore, we have 



V(x + x') = — — — sup E 

x (u,L)eA{x+x') 



\x + x' 



rrt I rrJ 

tXj iX/ _ _ , , 

< V(x). 

X 



and ( TTTT) gives the Lipschitz property of V. The proof of the Proposition is 
complete. 



4 Hamilton Jacobi Bellman Variational Inequality 



We have already seen that our value function V is increasing and Lipschitz 
continuous. These properties allow us to prove in a standard way that V 
satisfies the following Dynamic Programming Principle 

Principle 7 (DPP) 



V(x) 



sup E 

(u,L)eA(x) 



-r(tAr) 



V(XX' L ) + 



s dL< 



12 



for all t > 0, x > 0. 

For further literature on the subject, the reader to referred to Fleming, Soner 
(1993), Krylov (1980), or Yong, Zhou (1999) (theorem 4.3.3), for diffusion 
state processes or to Pham (1998) in the case of jump diffusion processes. 

We consider at this point the following HJB variational inequality: 

(14) 



m&x{H(x, V, V'(x)), 1 - V'(x)} = in R*_, 



V(0) = 0. 



where 



H(x,V,q) (15) 

= sup \ — rV(x) + axp(u)q + / [V(x — /(ax, y A u)) — V(x)] n(dy) \ . 
u>u I Jb J 



Let us recall that C 1,1 (1R + ) stands for the class of all real- valued, differentiable 
functions on IR + such that the derivative is locally Lipschitz. 

We also recall the definition of the viscosity supersolution, respectively viscos- 
ity subsolution. 

Definition 8 (i) Any lower semi- continuous (respectively upper semi-continuous) 
function v is a viscosity supersolution (subsolution) of ( fT^] ) if v(0) > (< 0) 
and 

max {H(x, ip, ip'(x)), 1 — f'(x)} < 0, 

(respectively > 0) whenever ip e C 1,:L (IR + ) is such that v — tp has a global 
minimum ( maximum ) at x > . 

(ii) A function v is a viscosity solution of if it is both super and subsolu- 
tion. 

Theorem 9 The value function V is a viscosity solution for the associated 
Hamilton- J acobi- Bellman Variational Inequality (fl^p . 



PROOF. First, we prove that V is a viscosity supersolution for (fl4"]l . In 
order to do this, let us consider x G W + and a C 1 * 1 test function ip such that 
V(x') — tp(x') > V(x) — (p(x) = 0, for all x' G Kl. Moreover, consider < h < x 
and the admissible strategy (u, L) G A(x) where L s = h, for all s > and u 
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is admissible and arbitrarily chosen. We have 



/ e~ rs dL s + e- r(thT) V{X 
J o 



(p( x ) = V(x) > E 
>h + E 

for all t > 0. We take the limit as t — > 0+ and get 

> /i + </?(x — /i). 

This latter inequality yields 

1 - ip\x) < 0. 



x,u,L\ 
tAr J 



(16) 



In order to prove H(x,<p,ip'(x)) < 0, we consider the admissible pair L s = 0, 
u s = u , for all s > (here Mo > w is arbitrarily chosen). We apply Ito's 
formula to e- r{t ^ V {X^ L ) to obtain 



E 



e-^V(X^ L )] - tp{x) 

tAr 



= E 
+ E 



rtAT . . 

/ (-re" r V W' ) + e- rs aX;' M ' L p(M )^(^'"' L )) ^ 

.J 

/ o MT jT e"" (<p (X X S ^ L - f (aX^, t/A«„))-^ (x^ L )) fi(dsdy) 



Recalling that (p(x) > E 



,- r {tAT) V (yx,uM 



V(X?£' ) , and dividing by t > 0, we have 



> E 
+ E 
> E 
+ E 

- O I E 



llT (- re ~ r VPC*"' L ) + e~ rs aX*> u > L p(u )<p'(X*> u > L )) ds 

- J y—rip(x) + e~ r *arEp(ito)v?'(:c)J ^ s 
1 r tA 



j_ rlAT r . . 

- y (e~ rt y? (a; - / (ax, y A w )) - V 5 ( x )) n{dy) 



sup e 

S<tAT 



X^ U ' L - X 



(17) 



where 0(5) ^ whenever <5 ^ 0. 



We wish to prove that E 
order to do this, we use 



sup s<tAr e 



j£x,u,L 



0, when t — > 0. In 



X x,u,L _ x 



< f ap(u )X*; u ' L ds' 
Jo 

+ f + I f(aX x s r> L ,yAu )^ds'dy). 

JO JB 
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Therefore, with the notation Cq = QdhM ^ we have, for some constant C, 



Co 

x x,u,L _ x 



< x (e Cos - l) 
+ Cx [ S+ [ e CoS 'fi{ds'dy) 

JO JB 



for all < s < t A r (we use the Lipschitz property of / in x uniformly in y, 
/(0, •) = and the upper bound for Xp u ' L given by (112]) ) . We multiply the 
last inequality by e~ rs , take the supremum over all < s < t A r, then the 
expectation with respect to P to obtain 



lim E 



sup e rs 

S<tAT 



X x,u,L _ x 



0. 



Notice that 

*!£M> 1 -p( T <t)>i-P{ m <t), 

where rji is the first time a claim occurs (it follows the exponential law). 
Consequently, 

lim = 1. (19) 

Returning to ( 1T7|) we let t — > 0+ and use (1T8|) and ( Tl9|) to get 

> {-rip(x) + axp^o)^'^)) (20) 
+ / if ( x ~ f ( ax i V A u o)) ~ V 0)} ir(<ijl) 

Combining ( 1201) and (116p . we prove that is a viscosity supersolution for (114p . 



In order to prove that the value function is a viscosity subsolution for (Tl4l) . 
we fix x > and consider an arbitrary test function (p G C 1 ' 1 such that 
y(a;') — (^(a;') < V(x) — (f(x) = 0, for all x' G R + . Let us suppose that the 
subsolution inequality does not hold. Therefore, there exists 5 > such that 

max {H(x, <p, (p'(x)), 1 — (p'(x)} < —5. 

We use the continuity of H and of cp' to obtain the existence of some rj G 
^0, x A j^-j , where K v denotes the Lipschitz constant for p on [0, e r x] , such 
that 

max {H(x', <p, p'(x')), 1 - p'(x')} < -5, if x' G B(x, rj). (21) 

Let us consider an arbitrary strategy (u, L) G A(x) and let X X,U,L denote the 
solution of (jHl) for (u, L) instead of (u, I). We define the stopping time 

(j = inf{t>0: X?' U ' L <£ B(x,rj)}. 
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Obviously a < r (the ruin time). We apply Ito's formula to e r ^ tAa ^ (p(X^ ,L ) 
and write 



E 





r ftAff , 


E 


IX 


[—re 












" rtAa r 


+ E 








Jo Jb 






■ rtAa 




E 


/ 






Jo 


+ E 


E e" 






S<tA<T 



re- rs <p{X*> u ' L ) + e- rs aX*> u ' L p{u )<p'{X*> u ' L ) ) r/.s 



(22) 



For s < t A a we have, from ([2T 



-re- r V(^' u ' ) + e- rs aX*' u ' L p(u )<p'(X*' u ' L ) (23) 
+e— jf (X?»> L - f (aX^ L , yAu ))-p (x^ L )) n(dy) < -fcT" , 



and, again from (J21 



It follows that 



<p' (X*> U ' L ) > 1. 



p (x:i"' L - AL S ) - cp (X X S ^ L ) < -AL S 



(24) 



Using (1251) we get 

£ jf (-re^VW' 1 *' 1 ') + e- rs aXf' ,1 ' L p(Mo)^ , (Xf '"' L )J 
+E [jT* jf e— ( v (x:l u ' l - f (aX^ L , yAu ))-p (x^ L )) »(dsdy) 
<5E 
-E 



-r(tAcr) ^ 



jT* j B e~ rs (<p (X^ L - f (aX^\ yAu ))-p (x^ L )) ix(dy)ds 



JO JB 

rtAcr r 

+ / ds e~ rs (p(x - f (ax,y Au )) - ip (x)) it (dy) 

JO JB 
rtho 



-E 
+E 

<5E 



JB 



— rs 



i.V (x ~ f (ax, y A u )) - p (x)) fj,(dsdy) 
£° J b e~ rs (<p (x:l u ' L - / (aX^ u ' L , yAu ))-p (x^ L )) fi(dsdy) 

1 1 — e ~r(tA<r)~ 

+ 4K v r]E . (25) 



3 -r{tAa) _ I 
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We return to ([22]) and use ([24]) and ([25]) to get 

~ e -r(tAa) _ Y 



E 



- r ^ V (XX' L ) - <p(x) < 5E 



+ AK v rjE 



I — e -r(tA<r) 



E 



tAo 



e~ rs dL, 



and, from this, 

V(x) = <p(x) 
> E 



^V(C L ) + / e~ rs dL 
Jo 



+ (6 - AK vV ) E 



1 _ e -r(tAa) 



(26) 



> E 



ip(X^) + e~ rs dL 



tAcr 



5 - 4K v r} 



E[t A a] 



for t small enough. We can suppose that a; is a strict global maximum point. 
Then there exists A > such that 

sup (V(x') - V (x')) = -A. 



We use ([2"B1) and write 

V(x) > E 



tAcr 



+XE 
> E 



e -r(tAv)y( X x,u^ + I e -rs dL 

e~ r{tA °h a<t 







+ s - iK *tP(* > t) 



r r t/Xa 
e-r^ViXX' ) + / e~ rs dL t 

+ (Ae-)A(^ f ) 



(27) 



The dynamic programming principle yields 



V(x) < sup E 

(u,L) 



e-^^yrxX' ) + / e~ rs dL 
Jo 



(28) 



Therefore, by the choice of rj < -J^- and A > 0, ( 1271) contradicts ([28]) . This 
proves that V is a viscosity subsolution for (|14p . Our Theorem is now complete. 



5 The Comparison Theorem 



The following Lemma provides an equivalent definition for the notions of vis- 
cosity super and subsolution. 
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Lemma 10 (i) A continuous function U is a viscosity supersolution for fl^D 
in if and only if, U(0) > and, for any x G IR+ and any test function 
ip G C 1 ' 1 sfic/i i/iai U — tp has a global strict minimum at x, we have 

max {H(x, U, <p'(x)), 1 - iff{x)} < 0. (29) 

(ii) A continuous function U is a viscosity subsolution for [Lfy in if and 
only if, U(0) < and, for any x G and any test function if G C 1 ' 1 such 
that U — if has a global strict maximum at x, we have 

max {H(x, U, <p'(x)\ 1 - <p'(x)} > 0. (30) 

PROOF. We only prove the assertion for viscosity supersolution, the proof 
for subsolution being similar. 

Suppose that (i) holds true. For any test function ip G C 1,1 such that U(x) = 
<p(x) and U — if has a global minimum at x, and all 5 > 0, we define 

(ps(x') = if {x) — 5\x' — x\ 2 , for all x G 

Then ip$ G C 1 ' 1 and U — tps has a global strict minimum at x. The assumption 
implies that 

max {H(x, U, <p' s (x)), 1 - ^(a;)} < 0. 

Obviously, U(x') — U(x) > (ps(x')-(ps(x), for all x' G R^_\{x} . The definition 
of H, together with the last inequality, yields 

max {H(x, ipg, ip'(x)), 1 — f'(x)} < 0. (31) 

Moreover, again from the definition of H, 

H(x,ip,<p'(x)) < H(x,ip s ,ip'(x)) 

+ sup ( / (ip (x - f(ax, y A u) - tp s (x - f(ax, y A u)) tv (dy) ) 

u>u \JB J 

< H(x,<p s ,<p'(x)) + C5, 

where C is a generic constant independent of S. We get, using (13T]) then taking 
the limit as 5 \ 0, 

max {H(x, ip, ip'(x)), 1 — f'(x)} < 0. 

For the converse, consider an arbitrary test function if G C 1 ' 1 and x G 
such that 

= U(x) - ip{x) < U(x') - ip{x'), 
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for all x' G R* + \ {x}. For e > such that e < f , we define 

5 £ = sup (U(x') - ip(x')) > 0. 

x'e-B(x,4e) 

It is obvious that lim £ ^ \ <5 e = 0. We introduce 

ip e = (U-ip- 5 £ ) l [0 , x -2e] + (17(0) - <p(0) - 5 £ ) 1 K _. 

We consider some sequence of mollifiers p n G (R; R+) , Swpp p n C B (o, 
and J R p n (t)dt = 1. Since C/ — </? is continuous, the sequence {p n * y? e } n con- 
verges uniformly on [0, x — 3s] to ip £ . Then there exists a subsequence (denoted 
by (p e )) such that Supp p £ C B (0, e) and 

[/ (a/) - v? (a;') - 25 £ < (p £ * <p £ ) (x') < U (x r ) - V {x') , 

for all < x' < x — 3e and all e > 0. Finally, we define the function 

F e (xf) = <p(xf) + (p e *<p e ) (xf). 

It is obvious that F £ G C 1,1 has the following properties: 

F £ (xf) = tp (x 1 ) , if x' > x — e, 
< U (xf) - 25 £ < F £ (x') , if < xf < x - 3e, 
F £ (x r ) < U (x') , if xf ^ x. 

The assumptions give 

ma3L{H(x,F e ,F' e (x)), 1 - F'(x)} < 0. 

Let us put 

G(x') = sup {-r (U(x') - F £ (x')) + ap(u)xf (<p'(xf) - F' £ (x')) 

u>u 

+ I (U(x' - f(ax', y A u)) — F £ (xf - f(ax', y A u))) n(dy) 
- I (U(x')-F £ (x'))n(dy)}, 

JB 

for any x' G R+. Then 

H(x, U, tp'(x)) - H(x, F £ , F' £ (x)) < G(x), (32) 

where 

G(x) < sup \ / (U (x - f(ax, y A u)) — F £ (x - f(ax, y A u))) ix(dy) 

u>u UB 

- l B (U(x)-F £ (x))n(dy)Y 
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We then consider the sets B u = jy G B : x — f(ax, y A u) G 5(x, 3e)| and get 

/ [U{x- f(ax, y Au)) — F £ (x - f{ax, y A u))) Tr(dy) 
Jb 

< f 28 £ 7r(dy) + Ctt(B u ) 
Jb\b u 

<2[38 £ + Cn(B u ), (33) 
where C > is a generic constant independent of e. Moreover, if y G B u , then 

x — f(ax, y Au) > x — 3e. 

Therefore, 

f(ax, y A u) < f(ax, y A u) < 3e. 

Since f(ax } y Au) > for xy > and f(ax, ■) is nondecreasing, we deduce the 
existence of some rj £ > such that r] £ —>■ as e — > and |/ G -B" only ii y < 1] £ . 
Thus, returning to (|33l) . we get 



(f/ (x - /(ax, y Au)) - F £ (x - /(ax, ?/ A u))) ir(dy) 
<C5 £ + Cn(Bn[0,r ]£ }). 

Consequently, 

G{x) < C5 £ + Cir{B n [0, %}). (34) 

Recall that ^ B. Thus, using (!3"4"]) in (132|1 and taking the limit as e — ■> 0, we 
obtain 

F(x, [/, < 0, 

and (i) follows. 



The assertion (ii) follows in the same way 

Under the assumption (A2) we are able to prove the following result on the 
comparison of viscosity solutions for (fT4"j) . 

Theorem 11 Let U and V be respectively a continuous viscosity subsolution 
and a continuous viscosity supers •olution for [Lfy both of at most linear growth. 
Then, if (A 2) holds true, we have 

U(x) < V(x), for all x G R* + . 



PROOF. For 5 > and e > 0, we denote by $ £j(S the function $ £i(5 : R + x 
K+ — ► R U {-oo} given by 

x ') = U{x) - V{x') - — (x - x'f -5{x 2 + (x'f) , (35) 
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for all x, x' > 0. Suppose that for some xq G K and some > we have 

U(x ) - 7(o;o) > 0. 

Since $ £j< 5 is upper semi-continuous and £7 and are of linear growth, there 
exists a global maximum point of $ e <5, denoted by (x £i s,x' eS ) G IR+ x R + . 
Obviously, since 3> £| <$(0, a/) < for all x' G 1R + , it holds that x £i s > 0. Moreover, 

Q 

7 £><5 = $ e ,5(x £>(5 , x' £tS ) > <S>e,s{ x o, x ) > 9 - 25x 2 > -, (36) 
for any 5 < S = Obviously, for 5 < 5 fixed, ( r y e ,s) E is increasing and 

If , \2 

j2e,S > le,S + ^ {X £ ,5 ~ X £ ,6 ) ■ 

Therefore, 

1 / \2 
Jl™- \^e,6-X efi ) =0. 

If, for all e > (or, at least for some arbitrary sequence e n such that e n — > 
when n — > oo) 5 = 0, then lim e ^ x e,s — 0, and, by taking the upper limit 
when e — > in (1331) . we get 

? < 17(0) - V(0) < 0, 

which contradicts the assumption 9 > 0. We deduce that, for s > small 
enough, and x' e s are strictly positive. We consider the test function 

ip{x) = v(x' EjS ) + y £ { x ~ X '^T + 5 ^ 2 + (<*) 2 )' for x g e;, 

such that U — if has a maximum point at x £t g. We write the variational in- 
equality and use the previous Lemma to get 

max {H(x E ,s, U, ip'(x £}5 )), 1 - if'(x £tS )} > 0. (37) 

In a similar way we have 

max {H(x' ej5 , V, ^«*)), 1 - ^« a )} < 0, (38) 

where 

ip{x') = U{x £)5 ) - ^( x e,5 ~ x ') 2 - S(xl s + {x') 2 ), for all x' G M* + . 

(a) We suppose that 

H(x £ ^ U, ip\x etS )) > H(x' £iS , V, ^ '«*))■ (39) 
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Then 



0< snv {-r(U(x £)S )-V{x' eS )) 

u>u v ' y 

+ ap(u) x e> s(p'{x et s) - x' s tp'(x' e s ) 



+ 



(U (x £i5 - f(ax EtS , y A u)) — V(x £S - f(ax' £S , y A u)fj ir(dy) 
+ J B (v(x' £;S )-U(x ej5 ))n(dy)}. (40) 

We use $(z £j5 ; x' e S ) > $(x £>(5 - f(x £iS , y A u), x' e S - f(x' £ S , y A u)) to get 
U (x £t5 - f(ax £t5 , y A u)) — V(x' £>s - f(ax' e S , y A u)) 

1 2 

U (x £tS ) - V(a£,*) + y ( Xe > 5 ~ /( ax ^' Z/ A m) — 4,i + y A M )) 



< 



2 £ x e ,S 



+ 5 [{x £) s - f(ax £7 s, y A u)f - x 2 eS + (x' £& - f{ax' eS , y Au)f - (< 5 ) 2 ) 
<U(x £>s )-V(x' £>s ), 

and, returning to (|4"U1) . we have 



< sup J -r ([/(x M ) - V(x' £ s )) + 2ap(u) 

u>u v ' 7 

< sup (2ap(w) - r) x (c/(x £)5 ) - 1/(4 5 ) 



2s 



+ 8 { xls + (x's 



Recall that sup u >„ 2ap{u) < ^giM < r ( see (A2)). Thus, it follows that 
7 £>< 5 < which contradicts fl36l) . 

(b) If ( BSD does not hold, we use fl37|) and fl38|) and we must have 

1-^(x £i5 )>0>1-^« 5 ), (41) 

thus 

We deduce that x' £ s = x £i $ = and get a contradiction. The proof of the 
comparison result is now complete. 



6 Numerical results 



We now turn our attention to some particular case and observe the optimal 
retention process by means of numerical simulation. We have seen that, for 



22 



the collective risk model introduced in [9], a single insurance contract is con- 
sidered and, of course, the risk is given for this one contract. A possible way 
to extend this model is to suppose that the risk concerns all contracts (or at 
least a percentage). We assume that the claims have constant intensity 6 and 
the random measure /i is associated with some Poisson process of constant in- 
tensity it(dy) = (3G${dy), where G$ corresponds to the Dirac mass. Moreover, 
the function / is given by f(x,y) = pxy, with < p < 1 (that is only some 
p part of the total contracts is subject to claims). In this case, the minimal 
retention level needed to cover expenditures is given explicitly by u = ( ' k2 ~ k ^ S 
and p{u) = (hi - k 2 ) (3p5+(l + k 2 ) f3pu, for all (fca ~ fcl)<5 < u < 5. 

Under the above assumptions, Eq. (jSj) reads 

X?' U ' L = x + a f X*> u ' L p{u s )ds - pa f X x s ^ L u s dN s - f dL s . (42) 
Jo Jo Jo 

Theorem 10 states that the maximized expected discounted dividends is the 
unique viscosity solution for the Hamilton- Jacobi-Bellman variational inequal- 
ity 

' max {H(x, V, V'(x)), 1 - V'(x)} = in R* , 
^ V{0) = 0, 
where 

H(x, V, q) = sup {— rV(x) + axp(u)q + j3 [V(x — apxu) — V(x)]} . 

(k 2 - kl )s <u<s 
k 2 - - 

The standard procedure in order to apply numerical arguments is to obtain 
a bounded space. Thus, we write the previous equation on [0, 1) by taking 
y = ^-j- and if)(y) = V(x). This leads to the following HJB equation 

' m a x{G(y^,^(y)),l-(l-y) 2 ^'(y)} = 0m [0,1), 
>(0) = 0, 

where 

G(y,iJ } q) = 



sup < —rip(y) + ap(u)y (1 — y) q + (3 

(k 2 -^)s <u<s { 
k 2 - - 



y(l-apu)-apu _ 
apuy + 1 — apu 



As in Mnif, Sulem (2005), the approximate solution of Eq. (|44|) is computed 
with the help of finite difference approximations and the policy iteration al- 
gorithm. 
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We consider two particular cases: the first one illustrates the natural frame- 
work in which the reinsurance company perceives a relative safety loading 
greater than that of the insurer, while the second example assumes the oppo- 
site. The data set we use is given in the following table 





h 


k 2 


5 


r 


P 


p 


Figl 


0.2 


0.25 


1 


0.07 


0.0011 


10% 


Fig 2 


0.2 


0.19 


1 


0.07 


0.0011 


10% 



For the first framework, the optimal retention level turns out to be maximal 
as shown by Fig 1. 



Optimal retention level for delta=1 (k1<k2) 



0.8 



o 

Cl 



2 0.6 



CD 



0.4 



Q. 

O 



0.2 



10 15 20 25 30 

Initial fortune 



Fig 1. Optimal retention level for 8 = 1, k% — 0.2, k 2 = 0.25 



As can be expected in the second case, if the initial reserve is great enough, 
then the direct insurer should play the safety card in order to maximize ex- 
pected discounted dividends. Indeed, since the relative safety loadings guar- 
antee a proportional steady income to the insurer, the optimal retention level 
is null (see Fig 2). 
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x iq- 16 Optimal retention level for delta=1 (k1>k2) 
4 1 1 1 1 

3.5 - 

3 - 

2.5 - 

2 f 

1.5 

1 

0.5 

o LI 

-0.5' 1 1 

50 100 150 

Initial fortune 

Fig 2. Optimal retention level for 5 = 1, k\ = 0.2, k 2 = 0.19 
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